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Abstract 

We present a detailed spectroscopic study of the hot gas toward the Galactic bulge along the 4U 
1820—303 sight line by a combination analysis of emission and absorption spectra. In addition to the 
absorption lines of Ovil K a , Ovil K^, Ovill K a and NelX K a by Chandra LTGS as shown by previous 
works, Suzaku detected clearly the emission lines of Ovil, Ovill, NelX and NeX from the vicinity. We 
used simplified plasma models with constant temperature and density. Evaluation of the background and 
foreground emission was performed carefully, including stellar X-ray contribution based on the recent X-ray 
observational results and stellar distribution simulator. If we assume that one plasma component exists in 
front of 4U1820— 303 and the other one at the back, the obtained temperatures are T = 1.7 ±0.2 x 10 6 K 
for the front-side plasma and T = 3. x 10 6 K for the back-side. This scheme is consistent with a hot 
and thick ISM disk as suggested by the extragalactic source observations and an X-ray bulge around the 
Galactic center. 
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1. Introduction 

X-ray surveys, such as the ROSAT All Sky Survey 
(RASS), show that there is a large enhancement of the 
emission in the Galactic center region expanding to ±30° 
in longitude and latitude, which is called an X-ray bulge. 
Snowden et al. (1997) summed 3/4 keV (0.4 - 1.2 keV) 
and 1.5 keV (0.7-2.0 keV) RASS data with 10° wide bins 
centered on £ = 353° and checked the latitude profile of the 
surface brightness. They found that the intensity distribu- 
tion of the enhancement in the b > 0° region and the b < 0° 
region significantly differed. They construct an isothermal 
(T = 10 6 ' 6 K) cylinder plasma model with an exponential 
fall-off in density with height above the plane and showed 
that the enhancement for b < 0° can be well explained by 
this model. The hot gas cylinder is located at the galactic 
center and its required radius is 5.6 kpc, with an electron 
density of 3.5 x 10 -3 cm -3 in the disk and a scale height 
of 1.9 kpc. They also suggest that in the b > 0° region, 
this model always produces less intensity than the data 
and this implies an additional component, such as Loop I. 
The total luminosity of the plasma cylinder was 1.9 x 10 39 
erg s _1 and the mass was ~ 3 x 1O T M . A molecular 
cloud of ~ 10 21 cm -2 hydrogen column density is opaque 
to soft X-rays lower than 0.5 keV. Shadowing observa- 
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tions using such molecular clouds located a few kpc away 
provide important clues. Park et al. (1997) observed the 
(£,b) = (-10°, 0°) direction with ROSAT, where a molec- 
ular cloud is located ~3 kpc away, and showed an anti- 
correlation of the molecular cloud density and X-ray in- 
tensity. They compare the X-ray intensity in on-cloud and 
off-cloud directions to find that almost 40% of the emis- 
sion in both the 3/4 and 1.5 keV band comes from behind 
the cloud. Almy et al. (2000) also carried out a shadowing 
analysis using ROSAT PSPC data for (£,b) ~ (337°, 4°), 
with an absorption column density assumed from IR (100 
fim) intensity, and revealed that more than 70 % of the 
emission is due to X-ray sources behind the cloud. 

Recent studies reveal that there are many kinds of X- 
ray emitters in our Galaxy. Not only the hot ISM in 
the Galactic disk and halo (Yao et al. 2009; Hagihara 
et al. 2010; Yoshino et al. 2009), but also CVs and nor- 
mal stars can contribute to the unresolved emission in the 
Galaxy (Masui et al. 2009). In order to distinguish be- 
tween the diffuse plasma emission and the sum of faint 
point sources, measurements of the absorption line due to 
the plasma is essential, and combination study of the ab- 
sorption and emission spectra can solve the density and 
scale of the plasma, as shown by the observation of the 
hot halo around our Galaxy using extragalactic objects 
as background sources (Yao et al. 2009; Hagihara et al. 
2010). 

Absorption observations toward the X-ray bulge were 
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carried out with an X-ray grating detector. The X-ray 
binary 4U 1820—303 is located in the globular cluster 
NGC 6624 at (£,b) ~ (2.8°, -7.9°), 7.6±0.4 kpc away 
(Kuulkers et al. 2003) as illustrated in figure 1. It was ob- 
served with the high resolution grating onboard Chandra 
and absorption lines of Ovil, OviII and NeVIII were 
reported (Futamoto et al. 2004; Yao & Wang 2006). 
Futamoto et al. (2004) analyzed data using curve of 
growth and constrained the column density of each ion 
as logTVovn = 16.2 - 16.7, logTVovm = 15.9 - 16.5 and 
logA^Neix = 15.7— 16.1. The velocity dispersion (v^) was 
derived using joint analysis of Ovil Ka and Ovil K/3 
to Vb > 200 km s _1 , and the column density of OviII 
and Ne IX are derived assuming the same velocity disper- 
sion constraint as O VII. Meanwhile Yao & Wang (2006) 
used their absline absorption model on the same data and 
reported the column density as logAovn = 16.3 ± 0.2, 
log N viii = 16.4 ± 0.2 and log7V Ne ix = 16.0 ±0.1 and the 
velocity dispersion was 255±J£ 4 km s- 1 . 

In this paper, we will use an emission spectrum obtained 
by Suzaku for a region ~ 2° away from 4U1820— 303. 
A combined analysis using the same method as used in 
the Galactic halo study with extragalactic objects by Yao 
et al. (2009) and Hagihara et al. (2010), in consideration 
with foreground diffuse and discrete emission will give us 
further information about the density and scale of the hot 
gas in the Galactic bulge region. 

^ 8.5 kpc > 



Sun GC 




4U 1820-303 



Fig. 1. Geometry of the sight line toward 4U 1820-303 and 
the Galactic plane 



2. Observation and Data Reduction 

The observation logs of the Chandra and Suzaku data 
are summarized in table 1. 4U1820— 303 was observed 
three times by Chandra, but we used only the HRC-S + 
LETG data to obtain the best energy resolution. The 
dataset and data reduction process is basically the same 
as in Yao & Wang (2006), but we followed the procedure 
presented in Yao et al. (2009) to extract the first order 
spectra of the LETG. 

Our Suzaku observations were taken with the CCD 
camera X-ray Imaging Spectrometer (XIS; Koyama et al. 
2007, Mitsuda et al. 2007 ). The XIS was set to the normal 
clocking mode and the spaced-raw charge injection (SCI) 
was applied to the data during the observations. We used 
processing data version 2.2.7.18 for the two observations. 
In this work, we used only the spectrum obtained with 
XIS1. 

We adopted the same data screening as in Hagihara 
et al. (2010), i.e., standard data screening in addition to 
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the exclusion of the thick atmospheric neutral oxygen col- 
umn density in the line of sight (Smith et al. 2007) and of 
the Solar wind charge exchange (SWCX) (Fujimoto et al. 
2007). We found that the counting rate was constant as 
a function of the neutral oxygen column density for the 
cleaned data. Thus, there is no significant neutral oxy- 
gen emission from the Earth's atmosphere in the filtered 
data. The solar wind intensity obtained with the Solar 
Wind Electron Proton and Alpha Monitor (SWEPAM) 
aboard the Advanced Composition Explorer (ACE ) was 
checked. We removed the time intervals when the proton 
flux in the solar wind exceeded 4 x 10 8 cm -2 s _1 (Masui 
et al. 2009). With this criterion, we discarded 19.7 ksec 
of the exposure and the net exposure time became 32.2 
ksec. 

There are no obvious discrete X-ray sources in the X- 
ray images for the 0.4-1.0 keV energy range. However, 
there would be many low luminosity sources in this di- 
rection, thus we applied a wavelet analysis using CI AO: 
wavdetect, to detect such point sources. We found 23 
point source candidates in the fields and removed a circu- 
lar region with 1' radius centered at the source position 
from the data. These candidates can be caused by sta- 
tistical fluctuations, and so we have to check them. We 
then compared the spectrum including the possible point 
sources and the spectrum excluding them and found there 
is no significant difference. To increase the photon counts 
and reduce the statistical error, we decided to use data 
including the point source candidates in further analysis. 

We constructed instrumental response files (rmfs) 
and effective area files (arf) by running the scripts 
xisrmfgen and xissimarfgen (Ishisaki et al. 2007). To 
take into account the stray light coming from outside 
of the CCD fov, we used a 20' radius flat field as 
the input emission in calculating the arf. We also 
included in the arf file the degradation of low en- 
ergy efficiency due to the contamination on the XIS 
optical blocking filter. The versions of calibration 
files used here were ae_xil_quanteff_20080504.fits, 
ae_xil_rmfparam_20080901.fits, 

ae_xil_makepi_20080825.fits and 
ae_xil_contami_20071224.fits. We estimated the non-X- 
ray-background from the night Earth database using the 
method described in Tawa et al. (2008). 

3. Analysis and Results 

3.1. Absorption spectrum analysis 

We first calculated the equivalent widths (EWs). We 
constructed a model with a power-law and narrow gaus- 
sians and fitted the spectra in narrow ranges as shown in 
table 2, and obtained EWs for each line using the eqwidth 
tool in Xspec. The errors in the equivalent widths are ob- 
tained as following; 1) we calculated the 90% error range 
of the normalization of the gaussian function. 2) Next, we 
calculated the maximum error of the EW (EW max ) us- 
ing normalization values of the best fit (A^est) and maxi- 
mum (N max ) as EW max EW/N hest x N max . Comparing 
the previous work, our obtained value for the Ovil Ka 
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Table 1. Observation logs 



Target name 


(*,&) 




ID 


Instrument 


Obs. Start date 


Exposure (ksec) 


4U1820-30 


(2.8°, 


-7.9°) 


98 


Chandra HRC-S+LETG 


2000-03-10 


15.1 


Bulge 3 


(1.3°, 


-7.5°) 


500001010 


Suzaku XIS 


2006-03-06 


51.9 (32.2)t 



t Value in parenthesis is the exposure time after data screening (see text) 



Table 2. Equivalent widths of absorption lines in 4U1820-303 spectrum 





VII Ka 


VIII Ka 


Ovil K/3 


Ne IX Ka 


Nex Ka 


Centroid (keV) 
EW (eV) 
EW (mA) 
Fitting Range(A) 


0.5734 

o-65±8:iS 

24.4±£* 
20-22 


0.6537 
0.69^ 

19.9^.8 
18-20 


0.6663 
1Q+ 016 

^ q+4.4 

18-20 


0.9223 

U.OZ_ 22 

12-14 


1.022 
<0.16 
< 1.9 
11.4-12.5 



(0-65^0.14 eV ) is smaller than those in Futamoto et al. 
(2004) (1.19±g;|5 eV) and Yao & Wang (2006) (1.06^ 
eV). We confirmed that this difference comes from the 
modeling of the continuum spectrum, with higher order 
diffraction photons (Yao & Wang 2006). We used only 
the 1st order spectrum of the HRC observation, and mod- 
eled the continuum spectrum in a narrow energy range, 
to obtain reliable values. 

To obtain the column density of ions and to estimate 
the total column density of the absorbing material, we fit 
all the lines simultaneously. We first fitted the contin- 
uum between 0.54 and 1.1 keV with a continuum model 
of a blackbody and a broken power-law. Foreground ab- 
sorption by the neutral ISM with solar metal abundances 
by Anders & Grevesse (1989) is assumed. This model 
returned a minimum ^ 2 /dof= 858.1/577 and the best-fit 
column density of the neutral ISM Nu is 1.9 x 10 21 cm -2 . 
This value is consistent with Nu estimated from by the 
100 /im intensity (2.3 x 10 21 cm -2 ) and reported in Yao & 
Wang (2006). 

This poor fit might be caused by residuals at the metal 
edge due to uncertainties in the Chandra response matri- 
ces. Nicastro et al. (2005) set the metal abundance of the 
neutral absorption material to be free to deal with this 
problem. As we followed their method, it decreases the 
residuals (x 2 /dof=761. 1/560) with the virtual absence of 
Ne and Na, and an over- abundance of C, N and Mg by 
factors of 3, 3 and 26 respectively. These values are due 
to the uncertainties in calibration, rather than true ISM 
metallicity (Nicastro et al. 2005). Though the x 2 is still 
not good, it is difficult to compensate the calibration un- 
certainties further and we use this model to describe the 
continuum. 

We next added two absorption lines representing O VII 
Ka and Ovil K/3 using the absem model, which is the 
same as in Hagihara et al. (2010). The column density 
and velocity dispersion of these lines are linked because 
both lines originate from the same ionization state of oxy- 
gen (model Al). The results are given in table 3. We next 
added absorption lines representing OviII Ka and NelX 
Ka step by step (model A2, A3). In each step, the ve- 
locity dispersion of lines are linked together. The results 



maintain consistency, thus in further analysis the velocity 
dispersion is always linked together. For the next step, 
we evaluated the hot phase hydrogen column density as- 
suming solar abundances of O and Ne (model A4) (see 
figure 2). We then let the Ne abundance vary and fit- 
ted the spectrum (model A4'). Please note that we fixed 
the abundance of O to the solar value, because the con- 
tribution of continuum emission is small and it is hard 
to determine the absolute metallically. Though the best 
fit hydrogen column density become smaller and the Ne 
abundance is 2.2ii"J, all the parameters are consistent 
with those of the A4 model. We will use this A4 model for 
further combined analysis, and evaluate the Ne/O abun- 
dance effects if necessary. There might be intrinsic ab- 
sorbing material around the 4U1820— 303. Futamoto et al. 
(2004) discuss the possibility of intrinsic absorption using 
a photo-ionization simulator and conclude that the size of 
the binary system and the luminosity cannot explain the 
O VII ionization fraction estimated by the column density 
ratio. Thus we consider that all of these absorption lines 
originate in the hot ISM in the further analysis. 

3.2. Emission spectrum analysis 

3.2.1. Modeling of the contribution from foreground 
emission and stars 

The emission spectrum from the 4U1820-303 vicinity 
obtained by Suzaku is shown in figure 3. The emission be- 
low 2 keV is very bright, and we can easily distinguish the 
emission lines of O VII, OviII, NelX, and NeX. We evalu- 
ated the intensity of these lines by a simple power-law and 
Gaussians fitting model (model El) and the results are 
summarized in table 4. Hereafter, the intensity of the line 
emission is shown in LU (Line Units), which corresponds 
to photons s _1 cm _2 str _1 . These line ratios are hard to 
reproduce with a single -temperature plasma. Thus we 
will study several components that can contribute to the 
emission in this region. We will consider the diffuse emis- 
sion from the Solar neighborhood, Loop I, and contribu- 
tion from stars in the Galaxy, step by step. 

As shown by the Suzaku shadowing observations of 
MBM-12 molecular clouds (Smith et al. 2007) and the 
evaluation of the soft X-ray background spectra of 14 fields 
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Table 3. Parameters of the hot ISM absorbing the 4U1820-303 spectrum 



Model 


Vb 




log 


N 




Abundance 


logT 


XVdof 




(km s _1 ) 


Ovu 


(cm 

viii 


- 2 ) 

Ne IX 


Hhot 


Ne/O 


(K) 
Ti 




Al 
A2 
A3 
A4 
A4' 


163 + ^ y 
239t?£ 

1Q1+254 
iyi -101 

112^ 4 
ISTtgg 4 


10.1/ _ .21 

i6.ioi° ;?^ 

16.141°;?* 


i5-9o±8:f? 

15.92lg;S 


15-87±g:i§ 


19.66t° ;{ 6 6 
19.46lg;S 


o 9+I.8 
z ' z -1.3 


u - z ' -0.07 
a. 90+0. 06 
°' Z ^-0.06 


708.4/573 
671.2/571 
759.9/569 
662.1/570 
659.8/569 



Table 4. Apparent surface brightness of each line, (model 
El) 



Ovu 


viii 


Ne IX 


Nex 


(LU) 


(LU) 


(LU) 


(LU) 


23.45iLn 


18.96t££ 


10.491S;S 


S QQ+0.76 

o.yy_o. 58 



(Yoshino et al. 2009), emission around our Solar neigh- 
borhood due to Solar wind charge exchange and/or the 
Local hot bubble is present. Most of the emission from 
the Solar neighborhood might be below 0.5 keV, except 
for ~ 2 LU of Ovil emission line. The apparent spectrum 
above 0.5 keV can be reproduced by an optically thin ther- 
mal plasma emission in collisional ionization equilibrium 
(CIE) with kT ~ 0.1 keV (Yoshino et al. 2009). As the 
uncertainty of the O VII line intensity is ~ 1 LU, we adopt 
these values as the foreground emission. 

The sight line toward 4U 1820—303 would go through 
the Loop I, a large structure seen in the radio waveband. 
This structure is considered to be an old SNR (10 6 years) 
located in Sco-Cen OB association. Egger & Aschenbach 
(1995) modeled this SNR based on RASS data. The den- 
sity and temperature in the cavity are 2.5 x 10 -3 cm -3 
and 4.6 x 10 6 K. We then assumed the column density 
and line intensity using these parameters as shown in ta- 
ble 5. After comparing these values with those in table 3 
and 4, we will neglect the effect of Loop-I in the absorp- 
tion line analysis, but take it into account in the emission 
line analysis. 

Recently, the contribution of stars in the soft X-ray 
background was studied using Chandra, XMM- Newton 
(Lopez-Santiago et al. 2007) and Suzaku (Masui et al. 
2009). The typical luminosity of a stellar corona is as 
small as 10 29 erg s _1 , but the number density is large in 
the Galactic plane and bulge direction. We estimated the 
contribution from stars by the following steps. First, we 
estimated the number of stars in the field-of-view, using 
a stellar population model. Secondly, we calculated the 
X-ray flux and spectra based on the currently available 
observational properties. 

We used TRILEGAL simulator ver 1.4 1 (Girardi et al. 
2005) to estimate the main sequence stellar distribution 
in the observing cone. This simulator synthesizes a stel- 
lar population for a given Galaxy field. Vanhollebeke 

1 http:/ /stev. oapd.inaf.it/cgi-bin/trilegal 



et al. (2009) used this simulator and compared the results 
with Two Micron All Sky Survey (2MASS) (Skrutskie 
et al. 2006) and Optical Gravitational Lensing Experiment 
(OGLE) (Udalski et al. 1997) observational data in 11 di- 
rections. There is a discrepancy of about 20% between 
the model and data, especially in the number of low lu- 
minosity stars. We used their best model parameters, to 
create the mock population, and counted the stars of sev- 
eral spectral types, and ages as shown in table 6. Clearly, 
there are lots of A, F, G, K, and M stars in the Suzaku 
field of view. Note that we count the stars in a radius of 
20' in the Suzaku sight line, and estimated the observed 
flux on the detector via the arf. Kuntz & Snowden (2001) 
estimated the X-ray luminosity of these stars based on the 
ROSAT data, and Rocks (2009) compiled the X-ray spec- 
tra after Lopez-Santiago et al. (2007) which uses XMM- 
Newton observation. We used the values shown in table 
7 as a template of each spectral type stars. In Lopez- 
Santiago et al. (2007), about half of the K and M stars 
are represented by the two-temperature coronal model, 
and the other half is described by a single -temperature 
coronal model. We added the flux from the mock star dis- 
tribution with these X-ray luminosity and spectra, by 8 
thin thermal plasmas as an empirical mock-up spectrum. 
The RS CVn type binaries are also bright in X-rays, with 
a typical X-ray luminosity > 10 30 erg s _1 . We adopt a 
simple exponential disk model by Ottmann & Schmitt 

(1992) , and found that the estimated numbers in a circle 
of 20' is about - 980. 

The HI column density toward the observing direction 
by 21 cm radio observation is 1.36 x 10 21 cm -2 (Kalberla 
et al. 2005). The HI +H 2 column derived by IRAS 100 
fim intensity with a conversion formula by Snowden et al. 

(1993) is 1.42 x 10 21 cm" 2 . We also calculated the hy- 
drogen column as a function of the distance from the Sun 
based on the global Galactic model by Ferriere (1998), 
and found that ^ 80% of the absorption material is lo- 
cated within 2 kpc from us. Thus, as a crude assump- 
tion, we do not apply absorption by the neutral ISM to 
the SWCX+LHB component, Loop I, and stars and RS 
CVns within 2 kpc, but we do apply the absorption of 
1.42 x 10 21 cm -2 for stars and RS CVns beyond 2 kpc. 
3.2.2. Hot ISM emission 

We then try to represent the observed energy spec- 
trum with contributions from stars and RS CVns and the 
Cosmic X-ray Background (CXB), which is modeled as 
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Table 5. Estimated contribution from the Solar neighborhood and Loop I 



5 



Model 




Line Intensity 






Column Density 








(LU) 






(log N/ (cm- 2 )) 






Ovil 


VIII Ne IX 


Nex 


Ovil 


VIII Ne IX 


Nex 


Solar neighborhood 


2 














Loop I 


1.83 


2.10 0.31 


0.04 


14.3 


14.8 14.3 


13.8 



Table 6. Criteria and number of stars in the direction of Suzaku observation within a radius of 20' for each spectral type 



spectral type 


(B-V)o 


My 


0-0.15 Gyr 


0.15-1 Gyr 


1-10 Gyr 


>10 Gyr 


O+B 


<-0.01 




8 


6 


7 





A 


-0.01-0.3 




6 


58 


12 


90 


F 


0.3-0.6 


2-8 


12 


78 


35212 


41152 


G 


0.6-0.8 


2-10 


6 


58 


97672 


29270 


K 


>0.8 


<8 


41 


240 


454739 


72546 


M (early) 


>0.8 


8-15 


133 


1054 


1851091 


407848 


M(late) 


>0.8 


>15 








35 






an absorbed power- law with a photon index of 1.4 and 
its normalization of about 10 photon sec -1 cm -2 str _1 
keV -1 at 1 keV (Hasinger et al. 1993). This model only 
reproduces about 1/10 of the emission below 1 keV. Even 
though we allowed the normalization of every stellar com- 
ponent to vary, it was impossible to exhibit enough O VII 
line because the temperature of the stellar components are 
all high. 

We next added a hot ISM component to reproduce the 
spectrum. We tried a one temperature hot ISM model 
with fixed and free (N, Ne and Fe) abundance and ob- 
tained poor x 2 / dof vame of 578.62/133 and 366.20/130 
respectively. This poor fit is caused by the spectrum 
where Ovil and NeX lines both exist and it is difficult 
for a single temperature plasma to reproduce these lines 
simultaneously. We then gave up the one temperature hot 
ISM model and tried a two temperature model with fixed 
abundance (model El). Though the x 2 was improved to 
358.17/131, there are still significant residuals between the 
model and spectrum as shown in figure 4. It is obvious 
that a continuum-like component is needed to reduce the 
residuals between 1 and 3 keV, thus we tried to free the 
CXB parameters. Though this model reduced the resid- 
uals and improved the x 2 /dof to 121.85/129, the best fit 
values of photon index is 2.4^q"i and the normalization 
is 41. photon sec -1 cm -2 str _1 keV" 1 at 1 keV. 
The flux is about 4 times larger than the nominal value 
(Hasinger et al. 1993) and is not reasonable for the CXB 
even though its fluctuation is taken into account. 

As mentioned before, there are large uncertainties in the 
number densities, luminosities , and the energy spectrum 
of the background low luminosity stars and binaries. For 
this reason, we tried the following two models. 

• Model E2: Assuming an underestimation of the 
stellar contribution, the normalization of the back- 
ground stars and RS CVns are set to be free. 

• Model E3: Assuming an unknown thermal compo- 



nent in the bulge, one high temperature thin thermal 
plasma emission model is added. 

In the E2 model, we set the normalization of stars to 
be free step by step and investigated the residuals and 
the normalization and estimated the effect of star contri- 
butions on the hot ISM properties. To compensate for 
the residuals around 2 keV, we needed an additional ~ 
2keV component. One plausible component to produce 
such hard emission is RS CVn binaries of which we mod- 
eled a spectrum with a two temperature plasma of 2.59 
keV and 0.17 keV. Moreover it is not well understood how 
many stars make binaries and there is uncertainty in the 
distribution of the RS CVn binaries. Thus we first set 
the normalization of background binary components to be 
free (Model E2-b: The emission model with background 
binary normalization is set to be free). Though we found 
that this model certainly reduced the residuals around 
2 keV, we also found that this model caused significant 
residuals around 1 keV and above 2.5 keV. Emission from 
plasma of kT =2.59 keV is a little too high to compensate 
for the residual. 

The M type stars exhibit coronal emission whose tem- 
perature is empirically known to be kT = 0.5 ~ 1.2 
keV, and characterized by complex emission with Ne, 
and Fe-L lines around 0.9 keV, (Lopez-Santiago et al. 
2007; Sciortino et al. 1999) and their low luminosity could 
cause large uncertainties in the number of stars within the 
FOV. For the next step, we fixed the binary normalization 
to the simulated value and set the normalization of the M 
type star to be free (Model E2-m: Emission model with 
background M type star normalization is set to be free). 
This model also could not explain the entire spectrum. 
However, residuals caused by these two models (E2-b and 
E2-m) are complementary. 

We thus tried to free the normalization of binary and M 
stars simultaneously (model E2-mb). This model fit the 
data with x 2 /dof of 144.37/131. The normalization ratios 
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Table 7. Stellar type parameters used in stellar emission estimation. The empirical spectral model are summarized by Rocks 
(2009) after Lopez-Santiago et al. (2007) and luminosities estimated by Kuntz & Snowden (2001). RS-CVn type binary parameters 
after Ottmann & Schmitt (1992) are also shown. 



Spectral 


Single 


Two temperature model 


Abundance 




log L x (erg s -1 ) 1 " 




type 


temperature 












Age 






kT (keV) 


kT hot (keV) 


kT cool (keV) 


EM cool / EM hot 


Z Q 


0-0.15 Gy 


t 0.15-1 Gyr 


1-10 Gyr 


F 


0.58 








0.5 


29.51 


29.02 


28.13 


G 


0.67 








0.3 


29.91 


28.82 


27.89 


K 


0.83 


1.17 


0.32 


1.97 


0.2 


29.66 


28.52 


27.89 


M (early) 


0.90 


0.80 


0.27 


2.02 


0.1 


29.31 


28.49 


27.54 


RS-CVn 


N/A 


2.59 


0.17 


0.22 


1.0 




30.75* 





f Luminosity in 0.1-2.4 keV 

I Luminosity in 0.4-4.0 keV independent of age 



to the simulated values are 6.57^2*14 an d 5.45+2.20 ^ or ^ 
stars and binaries respectively. It is difficult to confirm if 
these values are correct or not, so we accept this value at 
present. 

We then set the normalization of the K type star and 
F and G type star to be free sequentially (model E2- 
mkb, E2-mkfgb). This modification caused no significant 
changes to the fitting results and the contribution of K, 
F and G stars vanished. Finally we linked the normaliza- 
tion of all star components and set them to be free (model 
E2-(mkgf)b). The ratios of the obtained normalization to 
the simulated values are 2.79^0 78 an d 6. 22 +{'35 for all 
spectral type stars and binaries respectively. x 2 /dof is 
151.29/129 and it is hard to explain the whole spectra 
with this model. 

The normalizations and temperatures of the cooler and 
hotter hot ISM are consistent with each other in the three 
models (E2-mb, E2-mkb, E2-mkgfb). From these results, 
the normalization of the stars is of little effect on the hot 
ISM temperature. 

In the E3 model, an additional unknown component 
is assumed. We therefore fixed all the components to 
the simulated values, other than the hot ISM. First, we 
added a thin thermal plasma to the model because the 
contribution from an unknown stellar component would 
be highly possible (model E3). The abundance of the ad- 
ditional plasma is fixed to the solar value. However this 
model caused residual features like the E2-b model and 
could not fit the data (x 2 / dof = 165 - 37 / 131 )- The metal 
abundance of the stellar corona is not well understood 
and so we set the metal abundance to be free (model 
E3-A). This model explain the whole spectrum very well 
(x 2 /dof=121. 04/130). However, the best fit abundance is 
~ and the additional component is quite close to thermal 
bremsstrahlung of kT =1.1 keV. Actually we substituted 
thermal bremsstrahlung for the thin thermal plasma of 
E3-A model and found no changes in the fitting results 
(model E3-B). We then set the abundance of the plasma 
of E3-A model to 0.1 Solar value (E3-A'). This model also 
fit the data with x 2 / dof = 133 - 34 / 131 - 

Some 1.5 keV thermal components are reported in the 
spectra of nearby dM stars (e.g. van den Besselaar et al. 
(2003)) and this could be the origin of the unknown com- 
ponent. As shown in table 7, the metallicity of the coro- 



Table 10. Surface brightness of the O and Ne lines. 



Model OVII OVIII NelX NeX 

(LU) ( LU ) „ ( LU )„ ( LU ) 

E2-mb 14.72™ 9.80+^ 0.6lt^? < 0.19 

E3-A 7 15.80^;^ 10.75±^;|| 3.09^;^ 2.17^ 



nal spectra of the low luminosity stars is very low and 
this is consistent with the low metallicity of the unknown 
component. Stellar flares are another possibilitiy for the 
unknown component. In the flare period, the high tem- 
perature (> IkeV) component of the stellar corona be- 
come brighter than usual. It is not plausible that this un- 
known component originates from the hot ISM because 
the metallicity of the unknown component is not high 
(<~ 0.1) and the temperature and induced pressure is 
too high to maintain such plasma. 

From model E2 and model E3, we see that ~ 5 times 
normalization of stellar components or another thermal 
bremsstrahlung or a low abundance thin thermal plasma 
is needed to explain the whole spectrum. As shown in 
table 8 and 9, there is at most a 10% difference in hot 
ISM parameters between acceptable models. Our model 
assumes two extreme possibilities and it is reasonable to 
consider that the star contribution would change the hot 
ISM parameters at most by 10% in any case. Thus we use 
model E2-mb and E3-A' to represent the emission model 
in the further analysis. 

We evaluated the line intensities corresponding to the 
hot ISM component with model E2-mb and E3-A' as sum- 
marized in table 10. The intensity and ratio of the O lines 
are little affected by the stellar or additional background 
models, but the contribution of the Ne lines changes by 
the assumption of stellar components, especially of the 
high temperature stars. 

3.3. Combined analysis 

Two thermal components are at least required to de- 
scribe the emission spectrum. Assuming the geometry 
of the two temperature plasmas and the absorption to- 
ward the target, there are two possibilities for combined 
analysis. One is that only one plasma contributes the ab- 
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Table 8. Fitting results with model E2 (see text for details) 



7 



Model Foreground T ISM Background^" x 2 /dof 



cool hot Stars RS CVn CXB 
logT (K) Norm a logT (K) Norm a Type 6 Ratio c Ratio c 
[kT (keV)] [kT (keV)] 



E2-b 


fixed 


6.009+ u u « 
[0.088ir o i 5 ] 


215 


6.571™ 
[0.32llg;SS] 


26 


o + 2.Y 
"°-2.6 




fixed 


8 - 8i -1.28 


fixed 


167.78/132 


E2-m 


fixed 


6 019+ 067 
u.uia_ 074 

[0.090±°-°it] 


179 6 +212 - 9 


6 524+0- 035 

DlOZ4 -0.040 

i n no/i , 

[0-288t°'° 2 2 t] 


18, 


g+2.9 
'"-2.8 


M 


10.61±1;11 


fixed 


fixed 


163.12/132 


E2-mb 


fixed 


6-oi2±S-.gH 

[0.089±S:gl|] 


195.3lg^ 8 5 8 


6-5421°;^ 
[0.300t°;™] 


21 


.+3.0 
D -3.1 


M 


6.57l 2 ;JI 


5 45+ 1 - 27 
°- 4O -2.20 


fixed 


144.37/131 


E2-mkb 


fixed 


6-0121°;^ 

[o.o89±S:Si§] 




6.5431°;^ 
[0.30ltS:°&] 


21 


.+3.0 
°-3.0 


M 
K 


< 2.49 


4 93+ 1 - 48 

^• y<3 -1.43 


fixed 


142.54/130 


E2-mkgfb 


fixed 


6 016+ 060 
u.uiu_ 082 


183.4^ 


6.558lg;S; 
[0.312™] 


22, 


-1+3.4 
-2.8 


M 
K 
FG 


o 91 +1.67 

< 3.38 

< 1.23 


4 60 +1,63 

^• UU -1.35 


fixed 


140.42/129 


E2-(mkgf)b 


fixed 


6.oor±g;SI 

[0.088t°;™] 


204.2l 292 2 6 


6 519+ 027 
[0.2851°-™] 


21 


n+3.9 
U -3.6 


FGKM 


2 79+ - 74 

£ " ' -0.78 


fi 22 +1 ' 35 

°- ZZ -1.35 


fixed 


151.92/129 


E2-mb* 


fixed 


6 360 +0.020 

D.«3DU_ Q21 
|u.±tn _ Q oog j 


33.7^ 








M 


11 77+ 1 - 11 
±x - ' ' -1.00 


2 - 50 il.33 


fixed 


214.09/133 


E2-mkgfb* 


fixed 


6.334±S;ggi 


35.8+J-j 








M 
K 
FG 


7 ni +2 - 64 

'• U1 -2.17 

< 2.55 

4 29 +1-31 
^• zy -l.35 


4 04+1.52 

^•°^-1.59 


fixed 


196.99/131 



^All parameters are fixed to referred or simulation based value without mentioned 
^One hot ISM model. 

"Emission Measure 10 -3 J n e n p dl: in unit of cm -6 pc 



b Spectral type of the stars whose normalization is set to be free 
c Ratio of the normalization of stars of binaries to the simulation based value 



Table 9. Fitting results of model E3 



Foreground 
and backgrounds 



Model 



ISM 



Additional background 



X7ck>f 



cool 



hot 



logT (K) 
[fcT(keV)] 



Norm a 



+0.056 
0.062 
+0.012] 
-0.012J 
+ 0.078 
0.072 
+ 0.018] 
-0.014J 
+0.077 
0.071 
+ 0.018] 
-0.014J 
+ 0.061 
0.073 
+ 0.013] 
0.014J 



+ 182.1 



logT (K) 
[fcT(keV)] 

+0.023 
0.021 
+0.017] 
0.015J 

6.546l2;°£ 

[0.303tr o f 9 ] 
6 547+ - 026 

V-D41 _q 025 

[0.304+0.019, 



Norm a 



logT (K) 



Norm a 



Abundance 



+2.6 



+0.104 

0.082 
+0.816] 



+~27r~ 

2.0 



E3 
E3-A 
E3-B 
E3-A' 



E3-A'* 



fixed 
fixed 
fixed 
fixed 



6.019 
[0.090 
6.022 
[0.091 
6.023 
[0.091 
6.017 
[0.090 



205.1 



171 o +190 - 7 

1 ' x - u _88.4 



169.4: 



1QQ 7 +224.9 



6.567; 
[0.318 



27.3 



7.544 
[3.017 



17.4 



24.0 4 



23.9 



+2.8 



6.548 



-0.017J 

+ 0.026 



25.9 



+2.7 



109^ 



+0.015 
0.015 
+0.008] 



+2.3 
2.3 



[0.305 



+ 0.019] 



0.110 

109+0- 303 l 

• J - uy -0.248J 
^ + 0.108 
-0.098 
! +0.311] 
-0.222J 
7 ooq + 0.096 
' -^°^-0.077 
- + 0.372] 
-0.243J 



64.3 



+ 27.3 
16.5 



106: 

101"* 



69.6 



+ 26.2 
18.4 



38.5 



[1.495j 



+4.8 
5.8 



1.0 (fixed) 



<0.03 



0.1 (fixed) 



+0.038 
0.039 
991 +0.112] 
ZZ1 -0.104J 



d+4.7 
3 -4.7 



165.37/131 
121.04/130 
120.94/131 
133.34/131 



fixed 



6.440 
[0.237; 



37.3 



7.151 



0.1 



231.45/133 



^ All parameters of LHB+SWCX, Loop I, stars, and CXBs are fixed to referred or simulation based values 
*One hot ISM model. 

a Emission Measure 10 -3 J n e n p dl: in unit of cm -6 pc 
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Fig. 2. O and Ne absorption lines in the 4U 
1820—303 spectrum. The solid lines are model A4. 



sorption, and the other is that both plasma contribute to 
the absorption. Thus, we consider these two cases in this 
section. 

3.3.1. Uniform model with one absorbing plasma 

We assumed first as a combined model that only one 
plasma contributes to the absorption (model CI). The 
geometry is as follows; one uniform plasma (front-side 
plasma) exists in front of 4U 1820—303 and another uni- 
form plasma (back-side plasma) exists in back of 4U 
1820—303, as illustrated in figure 6. 

We constructed two isothermal plasma models with a 
uniform density, and a length L along the sight of line. We 




channel energy (keV) 

Fig. 3. Spectrum obtained in the galactic bulge region (left) 
and evaluation of the line intensities by Gaussians. (right) 

put an upper limit of 7.6 kpc on the length of the front- 
side plasma, to maintain consistency with the geometry. 
The velocity dispersions of the plasmas are linked together 
because this value was mainly determined using the ratio 
of Ovil Kq, to Ovil K/3 in the absorption spectrum. We 
tried four sets of models, Cl-1 and Cl-2 is the combination 
of E2-mb and A4 model, and Ne/O abundance ratio is set 
to be free in Cl-2. In Cl-3 and Cl-4 model, we adopt the 
E3-A' model for the emission model, and Ne/O ratio is 
set to be free in Cl-4. 

The fitting results are shown in table 11. The temper- 
ature of the front-side plasma (T~ 1.7 x 10 6 K) is deter- 
mined mainly by the absorption spectra and is consistent 
with the temperature obtained by the absorption analy- 
sis (T~ 1.8 x 10 6 K). The temperature of the front- and 
back-side plasma (1.7 x 10 6 K and 3.9 x 10 6 K) are both 
higher than the those determined only by emission anal- 
ysis (1.1 x 10 6 K and 3.6 x 10 6 K). Plasma of 1.7 x 10 6 K 
could emit Ovil and OviII lines three times and thirty 
times effectively than a plasma of 1.1 x 10 6 K. Though 
this makes the emission measure of the front-side plasma 
smaller to maintain the intensity of the O VII lines, O VIII 
lines are produced more effectively. Thus to suppress the 
O VIII intensity, the temperature of the back-side plasma 
becomes higher. Residuals caused by this adaptation 
could be compensated by the background components and 
this model also can reproduce the O, Fe and Ne emission 
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0.5 1 2 

channel energy (keV) 

Fig. 4. Spectrum and model (El). The dotted lines 
indicate foreground (unabsorbed) components and the 
solid lines indicate background (fully absorbed) compo- 
nents, hot ISM (green), CXB (black), SWCX+LHB, 
and Loop I (blue), stars except for M stars (or- 
ange), M stars (light blue), and RS CVns (magenta). 

lines. However in the energy range lower than 0.5 keV, 
there are residuals caused by the deficit of N lines, be- 
cause the lower temperature plasma (1.7 x 10 6 K) is too 
hot to emit N lines effectively. 

We summarized the physical properties of the plasma 
in table 12. We assume 1, 2 and 10 kpc for the length for 
the back-side plasma to estimate its density and pressure, 
because the back-side plasma contributes to the emis- 
sion, and we cannot determine the length and density 
separately. There could be two schemes. Assuming the 
back-side plasma confined in the Galactic bulge region, 
its length is at most 2 kpc and this leads to dense and 
high pressure plasma. Thus, this leads to a picture of 
that a hotter dense plasma existing around the Galactic 
center region and a warm thin plasma covering the disk. 
The other is that assuming pressure equilibrium between 
front-side plasma and back-side plasma, the length of the 
back-side plasma becomes ~ 8kpc, This means that a hot- 
ter plasma of large depth exists over the warm thin disk, 
because 4U 1820—303 lies 1 kpc below the Galactic disk. 
3.3.2. Uniform plasma model of two absorbing compo- 
nent 

For the second step, we assumed that two plasma com- 
ponents which exist in front of 4U 1820—303 contribute 
both the absorption and emission (model C2) as illus- 
trated in figure 7. Note that this assumes that there is 
no emission from the bulge behind 4U 1820—303. It is 
not determined whether the two plasma components are 
separate, or mixed together with some filling factors, but 
we again set an upper limit for the length to be 7.6 kpc 
and a common velocity dispersion. 

As with the previous model CI, we constructed two 
models with model sets of E2-mb + A4 (C2-1, C2-2) and 
E3-A' + A4 (C2-3, C2-4). The fitting results are shown 
in table 13. These values can be understood as follows; 
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Fig. 5. Best-fitted model and spectra by model E2-mb 
and E3-A' The dotted lines indicate foreground (un- 
absorbed) components and the solid lines indicate 
background (fully absorbed) components. hot ISM 
(green), CXB (black), SWCX+LHB, and Loop I 
(blue), stars except for M stars (orange), M stars 
(light blue), and RS CVns (magenta). Additional back- 
ground continuum in E3 model is indicated by gray line. 




Fig. 6. Schematic view of the model CI. 

First, two plasma components of temperatures ^ 3.5 x 
10 6 K and - 1.0 x 10 6 K are required to reproduce the 
emission spectra. The ionization fraction and emissivity 
of each plasma is determined by their temperature, then 
the column density and the emission measure to reproduce 
the absorption spectrum are obtained. We confirm this 
flow and found that the emission measures obtained here 
are about half of those obtained in the emission analysis. 
This is caused by the slight temperature decrease of hotter 
plasma induced by this combined analysis. 

With the model C2 results and the assumption of the 
Ne abundance of solar value, and the temperature depen- 
dence of the ionization fraction, the Ne X column density 
is at most 1 x 10 15 cm -2 . The upper limit of the column 
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Table 11. Fitting results of CI model 



Model 




Front-Side plasma 






Back-Side plasma 


Stars 


Binaries 


Additional 








at-,-. 


Length 


logT 


In e/ 




logT 


Norm a 


rtatio 


Ratio 


logT 


Norm 








(cm- 2 ) 


(kpc) 


(K) 




(km s _1 ) 


(K) 








(K) 






Cl-1 
Cl-2 


E2-mb 
A4 

E2-mb 
A4 


19 7 i+ u - uy 

t 

t 


4.67± a a ;™ 


fi 99 + 0.05 
D,ZZ -0.08 

t 

~ 99 + 0.05 
°- zz _0.06 

t 


1 4+ 1 - 2 

X - -0.8 
t 


109l^° 
126lf 5 4 


6.583r ; u r 9 

6.588± ;° ;i 




7.00±i;- 
6.67li;~ 








819.? 
819. ( 


Cl-3 
Cl-4 


E3-A' 
A4 

E3-A' 
A4 


19.68±°;^ 

19.53+-- 
t 


5-88^3.66 


6 20 +0 - 05 

°- ZU -0.04 

t 

R 90 + 0. 05 
°-^-0.05 

t 


t 


J -- L -34 

139+ 195 

— 55 


6 574+ - 041 

U - J '^-0.029 

6 592+ - 040 
o.oyz_ 04Q 


22.0^;° 
I8.9l«;i 






7 9M+ ' 085 
l .^_ 0g6 


37.7± 5 5 ; 2 

3T.it 5 5 ;° 5 


810.: 
805. < 



Note: All the parameters not written the table is fixed to the simulation based values. 
"Emission Measure 10 -3 J n e n p dl: in unit of cm -6 pc 
^Normalization of background M type star is set to be free 

c Ratio of the normalization of the background stars or binaries to the simulation based value 



Table 12. Physical properties obtained by the model Cl-4 



Model 


Component 


Length 
(kpc) 


Density 
(10- 3 cm- 3 ) 


Temperature 
(10 6 K) 


Pressure 
(10 3 cm- 3 K) 


Cl-4 


Front- Side plasma 
Back- Side plasma 


3.081?;- 
1 
2 

10 


fi+10.4 

4 s+ - 6 

^•°-0.4 
1+O.6 
°- i -0.4 


1 7+0.2^ 
L - 1 -0.2 
q+0.4 

°- y -0.3 

t 
t 


6.0±fe« 
16.8^-j 
12 1+ 3 - 7 

iZ - i -2.5 
k 1+1-3 
°- i -1.2 



but it effects little in the temperature of the ISM as shown 
by Hagihara et al. (2010). 

In table 15 we summarize the systematic errors of the 
background and foreground models and their effect on the 
best fit values of the hot ISM parameters. The obtained 
hot ISM parameters are all within 90 % statistical error 
and cause no change to our results and conclusions. 

In addition to the systematic errors mentioned above, 
we have to consider systematics caused by the combined 
analysis itself. Figure 8 shows contours of the tempera- 
tures of the hot ISM derived from emission analysis (E3- 
A') and combined analysis (Cl-4, C2-4). As shown in 
figure 8, contours only overlap partially and those from 
combined analysis are shifted to each direction. The best 
fitted x 2 /dof value of Cl-4 model is worse than those of 
the other models. These shifts would be caused by mod- 
eling a multi-temperature plasma with two temperature 
plasmas. To check the contribution of the uncertainty of 
the star and unknown component, we fixed the parameters 
of the unknown component to the best fit value of E3-A' 
model (table 9) and performed a combined analysis with 
the Cl-4 and C2-4 models. The results are consistent with 
those of Cl-4 and C2-4 model. These two models (Cl-4, 
C2-4) are representative of the other combined analysis 
models and we concluded the systematics originate in the 
combined analysis. Using these contours, the temperature 
range of the two plasmas are determined as 2. 8 — 4. 3x1 6 
K and 0.68- 1.9 x 10 6 K. 

We also checked the self-absorption effect of the emis- 
sion lines with the obtained column density and veloc- 



Sun GC 




Fig. 7. Schematic view of the model C2. 

density of Ne X from that of EW in the absorption spec- 
trum in table 2 is logA^NeX = 15.4, and is consistent with 
this upper limit. 

We summarized the induced physical properties of the 
plasma in table 14. The length of the hotter plasma is 
almost on the upper limit. While the cooler component 
has very short length < 0.7 kpc and the pressure is higher 
than that of the hotter by a factor of ~ 7. This assumption 
gives a scheme of a thin warm and a thick hot disk halo 
model. 

3.4- Uncertainty due to the model systematics 

In the absorption analysis, the energy resolution of the 
detectors (<~ 0.05 A) corresponds to ~ 400kms _1 and is 
not enough to determine the of the plasma only with the 
line shapes. The lower limit of can be determined from 
thermal limits. Thus, to determine the velocity dispersion, 
we used the ratio of the absorption depth of O VII K a and 
Ovil instead and assumed that the OviII, NelX and 
NeX originate from the same plasma, and linked the ve- 
locity dispersion of all lines. As and column density are 
coupled and if the plasma has a temperature gradient and 
local structure, this assumption causes systematic errors, 
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Table 13. Fitting results with C2 model 
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Model Uniform hot plasma 1 Uniform hot plasma 2 Stars Binaries Addition 







(cm- 2 ) 


Length 


logT 


Ne/O 


Vb 


^Hot 

(cm- 2 ) 


Length 


logT 


Ratio" 


Ratio 6 


logT 






(kpc) 


(K) 




(km s _1 ) 


(kpc) 


(K) 






(K) 


C2-1 
C2-2 


E2-mb 
A4 

E2-mb 
A4 


19.55ir 2 4 4 

19.55l° ° 2 4 2 
t 


76 +o.o 
76 +o.o 

' ,u -5.1 


6.48l U .r 4 
t 

fi zio + 0.04 
6 - 48 -0.04 

t 


9 +0 - 4 


139lf/ 
140tf 4 2 


t 

19 26+ 34 
t 


17+°-^ 

U - 1 ' -0.15 

n 1 s + 0.67 
°- 18 -0.16 


5 99 +u - 0V 

t 

5 99+ 07 

J - yy -0.15 
t 


7.58ii- 


4.73ii;iX 
4.56^;SJ 




C2-3 
C2-4 


E3-A' 
A4 

E3-A' 
A4 


19 57+ - 03 

-0.22 
t 

19 52+ 16 
t 


'■°-4.9 

5 9+ 1 - 7 


D - Oi _0.03 
t 

fi r.Q + 0.03 
°- DU -0.04 

t 


1 4+0-4 
X - -0.3 

t 


156l 2 6 f 
169l 273 


19 - 22 io.25 
t 

19.19±°;2 
t 


U,Z -0.1 
n io + 0.56 

u - id -o.n 


6 01+ - 07 

u,u -0.10 

t 

6 00+ 07 
t 






7 247+ - 069 

' -0.094 

7 97n + - 082 
' - z ' u _ .069 



Note: All the parameters not written the table is fixed to the simulation based values. 
a Normalization of background M type star is set to be free 

b Ratio of the normalization of the background stars or binaries to the simulation based value 
c Emission Measure 10 -3 / n e n p dl: in unit of cm -6 pc 



Table 14. Physical properties obtained by model C2-4 



Model 


Component 


Length 
(kpc) 


Density 
(10- 3 cm" 3 ) 


Temperature 
(10 6 K) 


Pressure 
(10 3 cm- 3 K) 


C2-4 


Plasma 1 
Plasma 2 


5.9+^ 
0-13±gif? 


1-8^ 

38.6± 4 3 f 5 3 


O+0.2 
°- z -0.3 

1 o +0 - 2 


38.6t§&° 



Table 15. Estimated systematic errors of the models and their effect on the best fit values of the hot ISM parameters (all the 
parameters are within statistical error). 



Component 


Systematic 


Affection to the best fit values of the Hot fSM 




Error 


Coof (T < 1.3 x 10 6 K) 


Hot (T > 2.5 x 10 6 K) 






Temperature 


Normafization 


Temperature 


Normalization 


Foreground stars 


+20% 


< 1%(+) 


-3% 


< 1%(-) 


< 1%(-) 




-20% 
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Fig. 8. 68%, 90%, and 99% confidence contours of the tem- 
perature of hot ISM for two parameters. Black, red and blue 
lines indicate the results of E3-A' model, Cl-4 model and 
C2-4 model respectively. The best fitted x 2 /dof values are 
indicated. 

ity dispersion of ions. When the column density O VII is 
10 16 cm -2 and the minimum velocity dispersion = 200 
km s _1 , the opacity r at the center of the resonance line 
is ~1.1, which is evaluated in the same way as Futamoto 
et al. (2004). The reduction due to the self-absorption for 
the resonance line is about 30% after integrating over the 
emission line profile. Note that we measured all resonance, 
inter-combination, and forbidden lines as O VII line, with 
the energy resolution of the Suzaku XIS. As the oscillator 
strength of the forbidden line is very small, the apparent 
reduction is diluted to be half as small as 15%. Also we 
neglected the effect of scattered-in photons from outside 
the line of sight. The r at the OviII line center is less 
than the unity with b = 200 km s _1 and A^ovm = 10 16 
cm -2 , and the total reduction is less than 20% without 
considering the O VII K/3 line. These effects will system- 
atically reduce the emission line intensity, but we neglect 
the effect in this paper because the correction is as much 
as the systematic and statistical errors. When we obtain 
the better energy resolution for the diffuse emission, we 
will able to evaluate by the ratio between the resonance 
and forbidden lines. 

4. Discussion 

In this paper, we analyzed emission and absorption 
spectra toward the Galactic bulge region, and found that 
there are at least two models to explain the emission and 
absorption data in this direction. We summarized the 
results in table 16. Only the emission measures are ob- 
tained for the back-side component, because there are no 
absorption data. As we simplified the models as much as 
possible, we will discuss the hot plasma within our Galaxy 
by possible extension from current simple models and by 
comparing the previous results at high latitude. 
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4.1. Extension of the current simplest geometry 

First, we took notice of the length of the front compo- 
nent. It is not likely that such a fully filled component is 
confined only in a sight line, and this implies that these 
components fill and prevail throughout the disk. Thus, we 
consider that such long extended plasma is a part of the 
hot ISM disk. Except for Plasma 2 in C2 model, the error 
range includes the upper limit determined by the geome- 
try, which is the distance of 7.6 kpc toward 4U1820— 303. 
It is not reasonable that the boundary of the two plasmas 
coincide at 4U1820— 303. It is useful to consider what 
happens when a plasma of the same temperature exists 
beyond 7.6 kpc as in figure 9. Additional plasma at the 
backside of 4U 1820—303, which does not contribute to 
the absorption, makes the emission measure of the orig- 
inal plasma decrease to maintain the emission intensity, 
while the column density is constant. The length becomes 
longer and the density becomes smaller assuming addi- 
tional distant plasma. If the density of absorbing plasma 
is reduced by a factor of 1/c, and the length is increased 
by the same factor c, the column density is kept as the 
same and the emission measure is reduced by 1/c. Thus 
if there is emitting plasma behind the target with den- 
sity n/c and length of c(c — 1) times the original one, the 
emission measure is recovered. It is easily written with 
the density n, the length L and a reduction factor c, 

CD = n/cxcL = nL (1) 

EM = (n/c) 2 x cL + (n/c) 2 x c(c - l)L = n 2 L (2) 

length=L, density=« 

• — ^ & 

length=cL , density=rc/c length=c( c- 1 )L, density -nic 

o — ^^^^^m&^^^m- 

Fig. 9. Schematic view of the extended plasma model be- 
yond 4U1820-303. 



Sun GC 




Fig. 10. Elongated models correspond to case CI and C2 
(Cl-B and C2-B) 

Both CI and C2 models can be extrapolated like figure 
10 as case Cl-B and C2-B. 



No. ] 

4-2. Comparison with the exponential disk observed at 
the high latitude 

The CI model implies a thin warm disk with a length 
of 3.1 kpc in the front and a hot plasma at the back, 
toward the bulge region. The hot plasma is possible to 
associate with the bulge, or to be a hot thick disk above 
the warm disk. We will compare the properties of the 
warm disk with the exponential disk model obtained by 
the combined analysis with absorption to extragalactic ob- 
jects; PKS2155-304 (Hagihara et al. 2010) and LMC X-3 
(Yao et al. 2009). In the extreme case, we can assume 
that the warm plasma has an extent of 7.6 kpc which 
corresponds to c = 7.6/3.1 = 2.5 in the above equations 
1 and 2. The total length and the total column density 
are calculated as (cL + c(c — 1)L) = c 2 L = 19.4 kpc and 
nL + (c — l)nL = cnL = 10 19 93 cm -2 . The thickness of 
the disk is 19.4xsin(6 4U ) = 2.7 kpc, where 6 4 u = -7.9° 
is the Galactic latitude of 4U 1820—303. This is compa- 
rable to that of the exponential disk. The column den- 
sity is also comparable to that found for PKS2155— 304 
considering the difference of the Galactic latitude, as 
10 i9.93 x s in(6 4 u)/sm(6pKs) = 10 19 " 17 cm" 2 , where b PKS = 
-52.2° is the Galactic latitude of PKS 2155-304. This 
value is consistent with the column density for the PKS 
2155-304 direction, log7V H = 10.10to.o7 obtained by an 
exponential disk model (Hagihara et al. 2010). Thus the 
front plasma of CI model corresponds to the exponential 
disk. 

To confirm this similarity, we applied an exponential 
disk model with parameters, which represent well both 
the emission and absorption spectra of PKS 2155—304 
(Hagihara et al. 2010). The exponential model well 
matched the absorption spectrum, but failed for the emis- 
sion spectrum due to the residuals below 1 keV. If we add 
components which only contribute to the emission like CI 
model, the fit returns a x 2 — 809.8/706 with two thermal 
component with logT = 6.5971^44 and 6 - 023 -o'.i83 with 
a stellar contribution the same as in E3-A'. 

The hot plasma at the backside was not observed at 
high latitude with LMC X-3 and PKS2155-304 (Yao et al. 
2009; Hagihara et al. 2010), but was detected in observa- 
tions toward other bulge regions (Almy et al. 2000). The 
temperature and estimated electron density assuming the 
size of 10 kpc is consistent with those by the RASS image 
model by Snowden et al. (1997). These support the idea 
that the back-side plasma is a hot plasma associated with 
the bulge region. The location, size, and pressure are, 
however, not determined by current observations, and the 
relation between the disk and the bulge is hard to be con- 
sidered. 

The C2 model implies a thin warm and thick hot disk 
halo model. Such a thin disk has not been observed by pre- 
vious studies, and its pressure is as high as 4 x 10 4 cm _3 K, 
and almost twice of the typical value at the midplane es- 
timated by Cox (2005) with thermal and non-thermal 
(Cosmic rays and magnetic field) components. The pres- 
sure of the warm plasma is also not balanced by a factor 
of 6. In addition, this simple C2 model assumes that the 
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hot ISM is spatially limited in front of 4U1820-303 and no 
plasma in the bulge region. If the hotter plasma can be 
extended beyond 4U1820-303 like C2-B in figure 10, the 
pressure will decrease in proportional to the density. We 
thus need some mechanism to confine the warm T ~ 10 6 
K plasma close to the Galactic disk. 

5. Conclusion 

We have analyzed high resolution X-ray absorp- 
tion/emission data observed by Chandra and Suzaku to 
determine the physical properties of the hot ISM toward 
the galactic bulge (4U 1820—303) direction with an esti- 
mate for the contribution from normal stars in the soft 
X-ray band. 

A two component plasma model can reproduce the ab- 
sorption and emission spectra. One model assumes that 
only one component contributes to the absorption in front 
of 4U 1820—303. The temperature, column density and 
length of the front plasma are determined as 1.7^q' 2 xlO 6 
K, 3.4l 2 ;^ x 10 19 cm" 2 and 3.1^ kpc. The temperature 
and emission measure of the back-side plasma are deter- 
mined as 3.9to'.3 xlO 6 K and 18.9t|£ cm" 6 pc. This 
model is consistent with a scheme with a hot X-ray bulge 
and an exponential disk model obtained from extragalac- 
tic source observations. If there are two plasma compo- 
nents contributing to the absorption, a thin warm plasma 
disk with temperature of I.O^q' 2 x 10 6 K and a length of 
0.13to.n kpc is confined to the Galactic disk, and it is not 
pressure-balanced with a hot thick disk with temperature 
of 3.2^; 2 x 10 6 K. 
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